Scaling and the prediction of energy spectra in decaying hydrodynamic turbulence 
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Few rigorous results are derived for fully developed turbulence. By applying the scaling properties 
of the Navier-Stokes equation we have derived a relation for the energy spectrum valid for unforced or 
decaying isotropic turbulence. We find the existence of a scaling function ip. The energy spectrum 
can at any time by a suitable rescaling be mapped onto this function. This indicates that the 
initial (primordial) energy spectrum is in principle retained in the energy spectrum observed at any 
later time, and the principle of permanence of large eddies is derived. The result can be seen as a 
restoration of the determinism of the Navier-Stokes equation in the mean. We compare our results 
with a windtunnel experiment and find good agreement. 



The Navier-Stokes equation (NSE) for hydrodynamic 
flow has been known for many years, and several interest- 
ing numerical results have been found. However, very few 
analytic statements have been derived from these equa- 
tions In this Letter we consider decaying isotropic 
turbulence, i.e. hydrodynamics without external forcing. 

One well known feature of the hydrodynamic equations 
is their invariance under re-scaling when the coordinates 
are scaled by an arbitrary quantity Any solution to 
the NSE can then be mapped onto another solution with 
a corresponding change in the velocity, the time and the 
diffusion coefficient v. The same scaling argument ap- 
plies to the energy density. In the following we shall 
show that this leads to a general scaling behavior of the 
energy density considered as a function of k,t,v. 

The result resembles the k 4 backscattering at integral 
scales obtained in EDQNM closure calculations S. In 
studies of the decay of homogeneous and isotropic tur- 
bulence self-similarity in the shape of the energy spectra 
are often assumed, expressed as the principle of perma- 
nence of large eddies (PLE) ^ . In a recent interesting 
phenomenological study it was proposed that for initial 
spectra not as steep as fc 4 there will be three ranges of 
self-similarity where PLE only applies to the very large 
scales [Q. 

We shall now derive the scaling properties of the en- 
ergy density directly from the scaling properties of the 
NSE. This has been done by one of the authors in the 
case where viscosity was ignored ||. Here we shall gen- 
eralize these results. Consider the energy density in D 
dimensions, given by 



£(k,t,u,K,L) = 

d D xd D y e ik(x " v 
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where L and K are infrared and ultraviolet cutoffs, re- 
spectively, and Q.d is the solid angle. We assume isotropy 
such that the energy density only depends on the modu- 
lus k — |k| of the wave vector. The total energy density 
is given by 



E(t) 



[ dkE{k,t,u,L,K) 

Jl/L 



d D x(v(xL,tf). 



(2) 



l/K 



Using the invariance of the unforced NSE under the scal- 
ing 



Ix, t->i 1_fc i, v^l h v, v^l 1+h v, 



(3) 



where h is an arbitrary parameter, the relation 
S(k/l, f-'H, l 1+h v, K/l, LI) = l 1+2h £(k, t, v, K, L) (4) 

trivially follows. 

In the following we assume the cutoffs such that for the 
relevant physical range we have 1/L •< k -C K. We will 
therefore suppress the explicit dependence of the energy 
density on the cutoffs. The viscosity provides a physical 
cutoff for large k while the cutoff 1/L is determined by 
the boundary conditions at the integral scale. 

Let us simplify (^) by introducing the new function ip 
by defining 



S(k,t,v) = k q ip(k,t,v). 



(5) 



In the above equations we introduced the quantity q — 
— 1 — 2h, which will be used in the following. To get more 
information, consider the relation Q as a function of I, 
differentiate with respect to I, and put I = 1 afterward. 
This yields the differential equation 



dh\k 
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= 0. 



(6) 



The general solution of this equation can be written in 
different equivalent forms: 

ip(k,t,u) = 

F 1 ( 2±3 in fc + In t, («+ 8 )+( 1 -g)T In fc + In t + 7 In v) (7) 

where F 7 is an arbitrary function, and 7 ^ is an arbi- 
trary constant. These solutions correspond to 
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£{k, t, v) = k q ip! (k^t, i/if^f ) 
for 7 = —(q + 3)/(l - g) or 

£(k,t,v) = k q ip 2 (k^t,vk 2 tj , 



(8) for large values of x. This gives 



(9) 



for 7=1. This is the main result of this Letter, the 
consequences of which we will describe in the following. 

If we take t = as the initial time, then the initial spec- 
trum is immediately obtained as a power k q as stated 
before (assuming that ip ^ for t = 0). If the cutoffs 
L, K had been explicitely included, they would appear 
as additional arguments kL and k/K in (||) and (||). 

If we simplify by assuming a very high Reynolds num- 
ber flow with a developed inertial range, then diffusion 
can be ignored for k and t sufficiently small such that the 
energy spectrum should have the form 



£ (k, t, v) 



k q V> ( k^t 



(10) 



Note that this corresponds to a solution (^) with 7 = 0. 
It is clear that in general the spectrum remembers the 
initial spectrum through the factor in front of ip, the 
scaling variable k^- 3+q ^ 2 t and the functional form of %p. 
The quantity k^ 3+q ^ 2 t is invariant under the rescaling 
(||) with the power h = — (1 + q)/2. We must stress here 
that ( |l0| ) is only valid as a 'local in spectral space' ap- 
proximation. Integrating ( |To| ) to obtain the total energy 
lead, as will be elaborated later, to an inconsistency. 

Shiromizu || has shown from renormalization group 
arguments that ( |l0| ) emerges after long time, due to the 
existence of a fixed point. 

From the scaling argument the actual form of the func- 
tion ip(x) cannot be found. If we impose the boundary 
condition that the Kolmogorov spectrum should appear 
for large fc's and/or t, we must have 



if)(x) 



, 3 + q 

k~t, 



(3D), 



for large values of x. This leads to 



K lo 9_5 + 3q_ 



(11) 
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For even larger values of x diffusion would of course be- 
come important. Now from classical K41 M dimensional 
counting we have £ ~ fc _5 / 3 e 2 / 3 , where e = dE/dt is the 
mean energy dissipation. Comparing with ( p"2| ) we have 
dE/dt ~ £ ( 5 + 3, ?)/( 3 +'?) which after integration gives 



E ~ f -2(l+q)/(3+g)_ 



(13) 



Note that for q < — 1 and ip(0) > we get from (|Io| ) an 
infrared divergency in the total energy. This reflects itself 
in an unphysical growth of the energy with time, so for 
q < — 1 ( |l3| ) must be modified by an explicit dependence 
on the cutoff 1/L. 

In 2D one expects the energy density to behave like 
k~ 3 , which corresponds to 



ip(x) 



(2D), 



(14) 



£ 



k-h- 2 



(15) 



for all q except q = —1, where £ = k~ 3 ip(t). 

The interpretation of the scaling as giving the devel- 
opment of an initial power spectrum is of limited interest 
in practice. In the following we shall therefore show that 
it is possible to use the scaling to obtain simple relations 
between the energy densities at different times, allowing 
any "initial" spectrum compatible with NS without forc- 
ing. Let us assume that the spectrum has been measured 
(or is otherwise known) at a time t ^ 0. From the scal- 
ing we then obtain 



£o(k,u) = £(k,t ,v) = k q ipi 



k 2 to 



(16) 



We do not assume any special form for the initial energy 
density £q. From the point of view of the scaling argu- 
ment, the quantity q is an arbitrary parameter. Later we 
shall explain how to determine q for a given spectrum. 

The relation ( |l6| ) can be used to construct the 
k— dependence of the scaling function from the measured 
£a by a simple shift of variables, 



(17) 



[k 2 t,vt 



k- q (i/i )-^5 £ (k (i/i )^,^). 

This in turn leads to the following relation between the 
energy densities at different times, 



£ k,t,v 



_2q_ 

3 + q 



£o k 



(18) 



This relation is a rewriting of the basic equation (4) . It 
is an expression of what could be called determinism in 
the mean in decaying turbulence, which indeed is not 
predictable in time. We see that in spite of the chaotic 
behavior of the velocity, it is in principle possible to pre- 
dict the energy density from knowing it at a fixed time. 
This predictive power is strongest in a region where diffu- 
sion can be ignored. Otherwise (18) makes predictions of 
the energy densities for different times and different vis- 
cosities, and relevant experiments may be hard to carry 
out in practice. 

On the right hand side of ( |l8| ) the time dependence of 
the argument proportional to v is absent for q = 1 and is 
very weak for q not too far from 1. In the application to 
be discussed below we see no effect of diffusion, except 
for large k. 

The scaling arguments leading to the above results are 
based on a scaling of the time by a positive quantity. 
Taking time itself to be positive, it follows that the ear- 
liest time we can have is t = 0. We shall now discuss 
how the parameter q can be determined in principle. In 
doing this we consider the high Reynolds number limit 
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by taking v — » (for q — 1 this is not necessary). Then 
( |l8| ) reduces to 

S(k,t,u^Q) = fa/t)*h £o(fc (i/i )^^^0) . (19) 

In the limit t — > existence of the spectrum on the left 
hand side requires 

£ (k (t/to)& ,v -> o) - (fc(t/to)*k)'. (20) 
In fll9] ) this implies 

£(fc,i,i/ -> 0) - fc 9 fort^O. (21) 

Thus we see that irrespective of whether the spectrum is 
measured at the initial time, it must behave like k q at the 
earliest possible time, in the limit v — ► 0, which means 
that k should not be too large. Even if the experiment 
was not operative at t = 0, one can extrapolate back to 
a hypothetical initial time t = 0. 
From ( po| ) it also follows that 



£(k,t,u -> 0) - k q for k -> 0. 



(22) 



Thus, for all times the spectrum keeps its original form 
for very low values of k. This is the principle of perma- 
nence of large eddies. It allows us to fix the parameter 
q if experimental results exist for sufficiently small k by 
fitting the data to k q . It should be noticed that we have 
only been able to prove the principle of permanency of 
large eddies when v — > 0, except for the case q = 1, 
where v is not mixed up with the time in ( |l8| ) . A weaker 
assumption used in the proof is the existence of the spec- 
trum for t — y 0, which imposes the boundary condition 
© on So. 

In order to test these results consider the experimen- 
tal measurements of the energy spectrum in decaying 
turbulence by Comte-Bellot and Corrsin ||. The en- 
ergy spectrum was measured in a windtunnel experi- 
ment at different distances from the turbulence gener- 
ating grid, see the figure (a) and (d)). The scaling pa- 
rameter q is determined from the data by using (13). 
It is found that E(t) ~ t~ XA which is in the range of 
most experimental findings ||. In panels (b) and (e) 
this is done and from the scaling we obtain in both cases 
q = 3.67. The panels (c) and (f) show the scaling func- 
tion iP(k^ +q ^ 2 t) = £(k,t, v)/k q obtained from (Q as a 
function of (fc(3+«V 2 i) 2 /( 3 +9). The collapse is seen to be 
perfect. The energy spectra show an inertial range of less 
than two decades, and even less for the 1 inch grid ex- 
periment. So even in this case of relatively low Reynolds 
number flow the dependence of the scaling function on v 
is weak and (|l^) holds very well. It should, furthermore, 
be noticed that the factor (i/t ) (1 ~ 9)/(3+<?) , which multi- 
plies v in Jl8| ) for these experiments vary from 1 to 0.42, 
i.e. by a factor less than 2.4. That could as well be the 
reason why the dependence on v can not be detected in 
the experimental data. 

In general, the energy cannot be obtained by integra- 
tion over the energy density, unless one has some knowl- 
edge of the scaling function ip. For example, without such 



knowledge one cannot judge the convergence of the rele- 
vant integral. In the case q = 1 the situation is, however, 
quite different, since (ph shows that 



£(k, t,u) = k ipi (k 2 t, v) = k ip (k 2 t) 



(23) 



where we left out any explicit reference to the diffusion 
coefficient since it is the same on both sides. In principle, 
this equation can be tested experimentally as a check of 
the validity of the NSE. 

Consider now the total energy, 



E(t) 



dk kip{k 2 t) 



1 



dy 1>(y), (24) 



which shows that in the case q = 1 the NSE predict a 
1 jt decay of the total energy. Since diffusion is included 
the integral is expected to converge rapidly. The inertial 
range, where the energy is constant, would be defined by 
introducing a cutoff yo(t) on the y— integral. Then one 
should require 



ya(t) 



dy ip(y) oc t 



(25) 



The Kolmogorov scale yo(t) can of course only be deter- 
mined from knowledge of the function ip. In order to ob- 
tain the energy in cases with q ^ 1, one would also need 
some knowledge of the scaling function. The classical 
problem of assuming a scaling relation of the form ( fioj ) 
where the diffusion is ignored all together is illuminated 
by considering the requirement of energy conservation. 
This leads to the exact relation 



dE 
~dt 



dk £{k,t) 



dt Jo 

= -2i> / dk k 2 £(k,t). 
Jo 



(26) 



Using (g|) after a shift of integration variable this can be 
expressed as 



2(1+9) 



q-1 

dx x 3 +i tpi( x i z ) 



, r°° „-i r 

/ dx x 3 +i zd z ipi(x, z) — 2z I 
9 Jo Jo 



(27) 

1 

dx xtpi(x, z), 



where the variable x is defined in (|1 
For q = 1 this simplifies to 



and z = v t 3 +<j 



dx ip\(x,v) = 2v j dx xtpi(x, v). (28) 
^0 

This relation gives a rather weak constraint on the scal- 
ing function. It is now seen why the second argument in 
■01 cannot be ignored: For q 7^ 1 (|27]) cannot be satis- 
fied, since this would lead to the meaningless requirement 
£(9-i)/(3+g) —constant. However, when the argument z 
is included there is no such problem, since all terms in 
( p7| ) are z— dependent. 

Eq. ( p8| ) simply states that (x) — 1/(2^). Thus, the 
mean value of the scaling variable is determined by the 
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diffusion, and (k 2 ) = \/(2vt). Thus, for iarge times, the 
mean vaiue moves toward small k, as expected, with a 
rate given by v. In a similar way we can analyze the 
content of (^7|). Defining the quantity 



((x a )) = dxx a i>x{x,z), 
Jo 



then (E7h has the solution: 



((jfe(«- 1 )/ 2 ))/((A ! ( 3 +9)/ 2 )) = const.z/t. 



(29) 



(30) 



This can be interpreted by introducing a weight function 
w(k) = fc^ -1 )/ 2 , from which we obtain 



(k 2 ) w ~ 1/vt, 



(31) 



so k approaches zero as time goes by, with a rate again 
given by v. It is interesting to note that ( |3l| ) behaves 
exactly as a spectral diffusion. 

To summarize, we have derived a rigorous scaling re- 
lation in the case of decaying hydrodynamic turbulence 
with a given initial state. This is expressed in terms of 
the function ip which itself is not a power function, but 
depends on the initial scaling. This means that the infor- 
mation of the primordial state is still on average present 
at later times. The information travels from the small 



scales to the large scales (or backwards in /c-space) in a 
diffusive way. 
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FIG. 1. Decay of grid generated turbulence in a windtunnel (Comte-Bellot and Corrsin, 1971)Panels (a) and (d) show 
energy spectra in two experiments with grid spacing of 1 inch and 2 inches respectively. The kinematic viscosity is that of air, 
v — 1.5 x 10 -5 m 2 /s. The spectra are measured at different distances from the grid, which by Taylors hypothesis is converted 
into times using t = D/Uo, where D is the distance from the mesh and Uo is the wind velocity upstream from the grid. The 
times are in panel (a): t=1.14 s (green), t=3.05 s (red), t=6.1 s (blue), t=9.8 s (purple), and in panel (d): t=1.07 s (green), 
t=2.49 s (red) and t=4.34 s (blue). Panels (b) and (e) show the total energy E(t) as a function of time t with the color coding 
as in panels (a) and (d). The straight lines are E ~ f _2 ( 1+ 9)/( 3 +9) w ;th q = 3.67. Panels (c) and (f) show the scaling function 
V>(fc (3+9)/2 t) = 8{k,t,v)/k q as a function of kt 2/{3+q) with q = 3.67. 
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